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Abstract 

We give a characterization of the central shell-focusing curvature singularity that can 
form in the spherical gravitational collapse of a bounded matter distribution obeying 
the dominant energy condition. This characterization is based on the limiting be¬ 
haviour of the mass function in the neighbourhood of the singularity. Depending on 
the rate of growth of the mass as a function of the area radius i?, the singularity may 
be either covered or naked. The singularity is naked if this growth rate is slower than 
R, covered if it is faster than R, and either naked or covered if the growth rate is same 
as R. 


1 Introduction 


An outstanding open problem in classical general relativity and relativistic astrophysics is 
to determine the hnal fate of the continual gravitational collapse of a star. The singularity 
theorems show that, given a few physically reasonable assumptions, such a collapse will 
terminate in a gravitational singularity. However, the theorems do not by themselves imply 
that the star ultimately becomes a black hole. In order for a black hole to form, the singularity 
must be invisible to far away observers. The theorems however allow for the possibility that 
the singularity may be naked, i.e. visible to a far away observer. Naked singularities, if they 
were to form in gravitational collapse, are expected to have theoretical and observational 
properties very different from those of black holes. Because of the difficulty in obtaining a 
general solution of Einstein equations, it is not known whether classical general relativity 
generically admits the formation of naked singularities in gravitational collapse |^. 
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Over the years, various examples of the formation of covered (i.e. not naked) and naked 
singularities have been found in general relativity, largely in studies of spherical gravitational 
collapse. Typically, it has been observed that for any given equation of state, both naked and 
covered singularities arise in these examples, depending on the choice of the initial density 
and velocity distribution of the star. While on the one hand there is a similarity amongst 
these various examples, suggesting an underlying pattern, on the other hand the issue of the 
genericity of these examples is an open one. 

In the present paper we give a characterization of covered and naked singularities in 
spherical collapse, on the basis of the behaviour of the mass function in the neighbourhood 
of the singularity. 

Our characterization is related to an elementary observation regarding light propagation 
in Newtonian gravity. The divergence of the density at the center of a collapsing spherical 
Newtonian star can be taken to mean that the ratio M/R^ diverges as i? —0, where M{R) 
is the mass to the interior of the radius R. However, according to the usual escape velocity 
argument, the escaping of light (moving at speed c) from this density singularity is governed 
by the ratio 2M{R)/R, not by the ratio M/R^. If in the approach to the singularity, 
2M{R)/R goes to a limiting value less than unity, light will be able to escape, but not 
otherwise. We may write 

2M _ 2M 

and it is evident that if M/R? diverges then 2M/R may either diverge or go to zero or be 
hnite and non-zero, depending on whether M/R^ diverges faster than 1/i?^, slower than 
1/i?^, or as fast as respectively. On the other hand, if 2M/R diverges, then M/R^ 

also diverges. 

Continuing with this Newtonian argument, let us consider the dependence of the mass 
function M{R) on the radius R, near the center, at any fixed time t. Clearly, M{R) will 
grow slower than i?°, and faster than or equal to R^, if we assume that the density at the 
center is non-zero, and if M{R = 0) = 0. If the evolution is non-singular, then M goes as 
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R^. A density singularity will form if there comes a stage during the evolution when M goes 
faster than R^. We know on the basis of theoretical and observational evidence that for a 
typical equation of state there are initial conditions (distribution of density and velocity) for 
which the collapsing matter cloud does not become singular, and either rebounces or attains 
equilibrium. These initial conditions correspond to selecting M behaving as R^ throughout 
the evolution. Let us now imagine choosing different classes of initial conditions, for the same 
equation of state, so as to move away from such non-singular evolutions and towards possible 
singular evolutions. In so far as M{R) is concerned, this is equivalent to changing the final 
behavior of M{R) from M going as R^ to M going as faster than R^, and approaching 
while possibly realizing growth rates between R^ and i?°. So long as we are considering the 
class of solutions for which the hnal growth rate is between R^ and R, the ratio 2M/R will 
go to zero and light will be able to escape the density singularity - this is an analog of the 
naked singularity in general relativity. If the growth rate lies between R and i?°, the ratio 
2M/R will go to inhnity and light will not be able to escape from the density singularity. 
This is an analog of the covered singularity. If M grows as R, either case may be realized, 
depending on the actual value the ratio takes. 

Again on the basis of theoretical and observational evidence, we expect that the black 
hole type (i.e. covered) singular solutions, which correspond to M going as faster than R, will 
be realized in gravitational collapse, from some initial conditions. As the initial conditions 
for a given equation of state are varied from the ‘weak’ to the ‘strong’, one goes from the 
non-singular class M ~ to the black hole type singular class (M grows faster than R). 
In doing so, one may or may not encounter solutions of the intermediate ‘naked singularity’ 
class (M goes slower than R but faster than R^). In this sense naked singularities, if they 
occur, he between non-singular solutions and covered singular solutions. 

Considering the similarity of spherical Newtonian collapse to the corresponding general 
relativistic case, one expects to be able to formulate such a characterization in spherical 
general relativity. We show below how this can be done, by considering the growth of 
a curvature invariant like the Kretschmann scalar, and by assuming the dominant energy 
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condition for the collapsing matter distribution. According to the dominant energy condition, 
pressure components cannot exceed the value of the density at any point in the spacetime. 

We also point out that we deal here only with the characterization of local nakedness of 
singularities, i.e. the issue of whether outgoing light rays at all terminate on the singularities 
in the past. We do not concern ourselves with the question of whether or not these light 
rays escape to inhnity (i.e. whether or not the singularities are globally naked). 

2 A classification of shell-focusing singularities 

We will assume the matter helds to have an energy-momentum tensor of Type I |^. With 
respect to an orthonormal basis [Eq, Ei, E 2 , E^), and with Eq being timelike, the energy- 
momentum tensor can be expressed in the form = diag{p,pi,p 2 ,P 3 )- All known matter 
helds, including those of zero rest mass, have an energy-momentum tensor of this form, with 
the exception of directed radiation, p is the energy density as measured by an observer whose 
world line has a unit tangent vector Eq, while the components Pa are the principal pressures 
in the three spacelike directions. The dominant energy condition implies that p > 0, |pq| < p. 

In comoving coordinates (t, r, 9, 0) the spherically symmetric line element in general rel¬ 
ativity is given by 

ds^ = E^dt^ — d^dr'^ — R^dVt^ (2) 

where cr, uj and R are functions of t and r. The components of the energy-momentum tensor of 
a collapsing matter distribution are given, in these coordinates, by = diag{p, —pr, —pr, —pr) 
The functions p, Pr and pt have the interpretation of being the density, radial pressure and 
tangential pressure of the matter held, respectively. 

If the matter content of the collapsing object is being described as a huid, then one must 
supplement the above description of the energy-momentum tensor with two equations of 
state, one each for the radial and the tangential pressure. If on the other hand the matter 
content is a fundamental held, say a massless scalar held, then the components p, Pr and pt 
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of Tl follow from the Lagrangian density of the scalar held. 
The Einstein held eqnations for this system are 


P = 


F' 

R^R'' 


( 3 ) 


R^R ’ 


(4) 


/ 2p; m' . 

= -+ ^^7 -^- -{PT - Pr), 


p + Pr Rip + Pr) 


(5) 


UJ = 


2p ARip + pt) 


p + Pr R{p + Pr) 


( 6 ) 


and 


e-B^ = I + /. (7) 

Here, prime and dot denote derivatives w.r.t. r and t, respectively. In eqnations (|^) and 
we have set SttG/c^ = 1. The fnnction F(t,r) resnlts from the integration of Einstein 
eqnations, and has the interpretation of being twice the mass to the interior of the shell 
labeled by the coordinate r. The fnnction /(t, r) is dehned by the relation e‘^ = R'"^ /(1 + /) 
and satishes the condition / > — 1. 

We assnme that at time f = 0 the spherical object starts nndergoing gravitational col¬ 
lapse, and we choose the initial scaling i?(0,r) = r. Next, we assume that a shell-focusing 
curvature singularity, given by R{ts{r),r) = 0, forms during the evolution. At this sin¬ 
gularity, the Kretschmann scalar Rabcd.R^'^'^'^ diverges. The shell labeled by the comoving 
coordinate r becomes singular at the time fs(r). In particular, the ‘central’ curvature sin¬ 
gularity, i.e. the one at r = 0, forms at the time ts(0). It is the central singularity that is 
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of primary interest from the point of view of nakedness, and is the object of study in the 
present paper. We will assume that fs(0) < ts{r ^ 0) and that the central shell-focusing 
singularity is not preceded by any shell-crossing singularity. Since F{t,r) is twice the mass, 
it is physically reasonable to take F{t, 0) = 0, prior to singularity formation. Furthermore, 
since the center is the first point to become singular, according to time t, we assume that 
F(f,(0),0) = 0. 

Let us consider the propagation of an outgoing radial null geodesic congruence C- As 
is known, the occurrence of a covered or locally naked singularity is determined by the 
behaviour of the geodesic expansion 6 = in the approach to the singularity. If the 
expansion of the outgoing null congruence continues to remain non-negative in the approach 
to the singularity, the singularity will be naked. On the other hand, the singularity will be 
covered if 9 becomes negative in the approach to the singularity. It can be shown [Q that in 
a spherically symmetric spacetime 9 can be written, using the form (Q) of the metric, as 

Here, R' and C,'^ are non-negative quantities. We will assume that the ratio 2R'(^/R 
remains non-zero in the limit. While this is an untested assumption, it is seen to hold at 
least in those cases when the dependence of R{t, r) on r is a power law - since that gives 
R'/R ~ 1/r, and also, = dr/dk ~ r/k. It then follows that the necessary and sufficient 
condition for the singularity to be naked is that F/R approach a value less than or equal 
to unity along the outgoing geodesic as it meets the singularity in the past. (The only 
exceptions to this condition arise when / takes the value —1 or -|-cx) in the limit. In the 
former case, the singularity could be covered even if F/R goes to unity in the limit, and in 
the latter case it could be naked even if F/R takes a limiting value greater than unity.) 

The above condition on F/R is nothing but the usual condition of no trapping, and is 
analogous to the condition on 2M/R in the Newtonian case. 

Consider next the rate of divergence of the Kretschmann scalar at the central curvature 
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singularity, as the singularity is approached in the past along an outgoing null geodesic. 
Because of the dominant energy condition, the leading divergence in the Kretschmann scalar 
will be due to the energy density p. Also, if the density is non-singular, this scalar will be 
hnite. 

We hrst transform from the comoving coordinates (t,r, 6*,0) to coordinates {t, R,6,(j)), 
thus using the area radius as one of the coordinates. From equation (|^) it follows that we 
can write, in these new coordinates 


p{t, R) 


{dF/dR)t 

R^ 


(9) 


In order to evaluate the growth of p along an outgoing null geodesic, we use the trajectory 
of the outgoing null geodesic itself as one of the coordinates. We label the trajectory by a 
parameter X, which will be a function of t and r. We eliminate the coordinate t in favour 
of X, and in the (X, R, 9, 0) coordinates we get the evolution of the density to be 


,^ {dF/dR)x + {dX/dR)t{dF/dX)n 
P[X,R) = -- 


( 10 ) 


Now, if the density evolution is to be non-singular, F{X, R) must grow as R^ or slower, 
as one approaches R = 0, keeping X hxed. (A faster approach rate will cause the density to 
diverge, as is evident from (llO|)). On the other hand, if the evolution results in a singularity, 
we cannot conclude anything about the growth rate of F as a function of R, at a hxed X. 
This is because the dominant divergence may come from the second term on the right hand 
side in Eqn. (pT]), rather than the hrst term. However, we know from that in order 
for the singularity to be covered, it is necessary that this growth rate be at least as fast as 
F(X, R) ~ R, or faster. 

We can now give a classihcation of the nature of the evolution (i.e. whether it is non¬ 
singular, or if singular whether it is naked or covered), depending on the behaviour of the 
mass function F{X, R) as a function of R, for a hxed X, in the approach to the singularity. 
For a given form of matter, (say a huid with a known equation of state, or a massless scalar 
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field) first consider initial data for which the evolntion is non-singnlar. This means that 
along an ontgoing nnll geodesic, F grows as or slower. 

Consider next changing the initial data so that we come to the class of data that lead 
to singnlar evolntion. We divide snch data into two possible snbclasses, one for which the 
limiting growth of F, for a hxed X, is or slower, and the other for which it is faster than 
R^. The singnlarities in the former class wonld be naked, since the ratio F/R will go to 
zero in the limit. Hence, if the growth of F{R) is as R^ or slower, the evolntion is either 
non-singnlar or singnlar and naked. 

Those singnlar evolntions for which the limiting growth rate of F is faster than R^ can 
be fnrther snbdivided into three classes. If F grows faster than R^ bnt slower than R, the 
singnlarity is naked. If F grows as R the singnlarity is covered if the limiting valne of F/R 
exceeds nnity, and naked if this limiting valne is less than or eqnal to one. If F grows faster 
than R the singnlarity will be covered. 

In view of the above characterization one can nnderstand better the conditions that 
are reqnired for naked singnlarity formation in spherical inhomogeneons dnst collapse. For 
simplicity, we recall the case of marginally bonnd dnst collapse, for which the initial density 
prohle near the center, written in a series as 

p{R) = Po + piR + + ••• (11) 

2 D 

determines the ontcome of the collapse. It tnrns ont that the singnlarity is naked if pi < 0 or 
if Pi = 0, p 2 < 0. It is covered if pi = p 2 = ps = 0. If pi = p 2 = 0 and ps < 0 the singnlarity 
is naked if the dimensionless qnantity ( = i/Sps/dpo^^ is less than or eqnal to —25.9904, and 
covered when ( exceeds this valne 0. 

One conld ask why the transition from the naked to the covered case takes place at the 
level of the third derivative. A physical answer is the following. It can be shown that in the 
approach to the central dust singularity along an outgoing null geodesic, R ~ r“, with a 
taking the value (1 -|- 2p/3) where p is dehned such that the hrst non-vanishing derivative in 
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the above expansion for the density is the pth one. According to our naturalness argument 
above, if the singularity has to be naked, F/R^ must not vary faster than Noting that 

F ~ r^, we can write F/R^ /^3(i-a)/a implies that a < 3, i.e. p < 3, is necessary for 

nakedness (p < 3 is sufficient for nakedness). This physically explains the signihcance of the 
transition taking place at the level of the third derivative. We also note that the condition 
that 2R'C,^/R does not vanish in the limit is satished in the dust model. 

3 Discussion and Conclusions 

We summarize here the assumptions that we have made: (i) We have considered Type I 
matter helds obeying the dominant energy condition; (ii) we have assumed that the chosen 
initial conditions are such that any central shell-focusing singularity that might form in 
collapse is not preceded by a shell-crossing singularity; (iii) in the expression for 9 we 
have assumed that the quantity 2R'(^/R remains non-zero in the limit. 

Subject to these assumptions, and subject to choosing a specihc form of matter, we have 
given a partial classihcation of the nature of evolution (as to whether it is non-singular, and 
if singular whether it is covered or naked) in terms of the behaviour of the mass function in 
the neighbourhood of the origin. For a given set of initial conditions, the evolution is non¬ 
singular or naked singular if this growth rate is R^ or slower. If at some epoch during the 
evolution the growth rate becomes faster than R^, this shows the formation of a curvature 
singularity at that epoch. This singularity is naked if the growth rate is slower than R, it 
is covered if the growth rate is faster than R, and it could be either covered or naked if the 
growth rate is same as R. 

We make the following observation about the nature of the shell-focusing singularity. 
Suppose it is the case that for a given form of matter, there are generic initial data leading to 
non-singular evolution, and other generic initial data leading to covered singularities, but the 
initial data leading to naked singularities is non-generic. This implies the following regarding 
the behaviour of the mass function. There are generic initial conditions for which the mass 
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function grows as and other generic initial conditions for which, near the singularity, it 
grows as R or faster. However, the initial conditions for which it grows faster than R^ but 
slower than R near the singularity are non-generic. 

A priori, one cannot draw conclusions as to how the dependence of the mass function 
on R changes, when the initial data is changed so as to go from the non-singular class of 
solutions to the singular class. If this change is a continuous one, then there will necessarily be 
generic data for which M{R) grows faster then R^ but slower than R, and naked singularity 
formation will be generic. On the other hand, if the change is a discontinuous one, and there 
is a jump from M{R) ~ R^ to M{R) ~ faster than R, naked singularities will be non-generic 
(provided also that evolutions of the type M ^ R^ are generically non-singular). Whether 
or not this change is a continuous one cannot be decided from the present study, and further 
investigation is necessary. 

Although we cannot conclude here about the genericity or non-genericity of locally naked 
singularities, we compare our results briefly with other recent developments relating to naked 
singularities in spherical gravitational collapse. Christodoulou [Q has shown that for the 
spherical collapse of a massless scalar held, the initial data leading to globally naked singu¬ 
larities is non-generic. To our understanding, Christodoulou’s result does not have a bearing 
on whether or not the locally naked singularities in his model are non-generic. Since the 
arguments we have given above are only for local nakedness, they cannot be directly com¬ 
pared with Christodoulou’s conclusions. Joshi and Dwivedi have demonstrated that the 
formation of locally naked singularities is generic in the spherical collapse of Type I matter 
helds obeying the weak energy condition. Since the weak energy condition does not nec¬ 
essarily imply the dominant energy condition, their result does not necessarily imply that 
M{R) grows faster than R^ and slower than R near the singularity, for generic initial data. 

Finally, we comment briehy on numerical studies of spherical collapse. In these studies 
(including those on scalar helds and huids) it has been found that dispersive evolutions and 
black hole formation both result from generic initial data. In many of the the black hole class 
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of solutions, the amount of mass that collapses to form the black hole behaves as {p — 
where p is a parameter labeling a family of solutions, and p* is the critical value below which 
the evolution is dispersive. The limiting solution for which p —> p* is a naked singularity, and 
numerical studies suggest that such a naked singularity is non-generic, because it appears to 
occur only at one value of p, i.e. p = p*. 

To our understanding however, the following important issue should be taken note of. In 
most numerical studies to date, the identihcation of a black hole in a simulation is carried 
out via the detection of an apparent horizon. Such a detection does not by itself rule out 
the development of a Cauchy horizon and a naked singularity which is not covered by the 
apparent horizon. Hence, in principle, the solutions that are classihed as black holes could 
include, as a subset, naked singular solutions. The few simulations that have been carried out 
using double null coordinates that can be extended up to the center of spherical symmetry 
do not appear to offer a completely clear picture in this regard. Hence it appears to us that 
further numerical investigation is necessary in order to examine the issue of genericity of 
naked singularities in spherical collapse. 
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nologica (JNICT) Portugal, under contract number CERN/S/FAE/1172/97. 
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